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Abstract. The two-kaon interferometry at RHIC is studied in a multi-phase 
transport model. Similar to the pion case, we find strong space-time correlation 
at freeze-out for the kaon emission source, which results in a large positive R ou t— t 
term and tends to reduce the Rout/R s idc ratio. Unlike the pion case, the source 
radii for kaons determined from the emission function are close to the radius 
parameters extracted from a Gaussian fit to the correlation function, ft- 



1. Introduction and Summary 

Particle interferometry based on the Hanbury-Brown Twiss (HBT) effect has been 
used extensively in heavy ion collisions to extract the information on the emission 
source of particles |2 El E]- In particular, the long emission time as a result of 
the phase transition from the quark-gluon plasma to hadronic matter in relativistic 
heavy ion collisions may lead to an emission source which has a much larger radius in 
the direction of the total transverse momentum of detected two particles (i? ut) than 
that perpendicular to both this direction and the beam direction (i? s ide) Since 
the quark-gluon plasma is expected to be formed in heavy ion collisions at RHIC, it 
is surprising to find that the extracted ratio i?out/-Rsidc from a Gaussian fit to the 
measured two-pion correlation function in Au+Au collisions at ^/s = 130A GeV is 
close to one [E1II1IEJ, very different from predictions of hydrodynamical models @JE]. 

Since particle interferometry probes the phase-space distributions of particles 
at freeze-out, it is natural to apply transport models to HBT. The reason is that 
particle freeze-out is dynamically generated in transport models when the mean-free- 
path exceeds the system size at later stage of expansion, whereas freeze-out has to be 
imposed in hydrodynamical models. Using a multi-phase transport (AMPT) model, 
we have found that the small pion i? ut / -Rside ratio could be due to a large and positive 
space-time correlation in the emission source [0]. Furthermore, the pion source at 
freeze-out is highly non-Gaussian, leading to much larger pion source radii than the 
radius parameters from a Gaussian fit to the three-dimensional correlation function. 

In this study we extend the work of Ref. |Hj by studying the kaon interferometry 
in central Au+Au collisions at RHIC energies. Using the AMPT model, we find 
that, unlike the pion case, the kaon source radii extracted directly from the emission 
function are close to the fitted radius parameters extracted from a Gaussian fit to the 
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three-dimensional correlation function. Our results also show that the kaon emission 
source has a large and positive correlation between time and position along the out- 
direction at freeze-out, similar to what we have found earlier for the pion emission 
source. We expect that the study of kaon interferometry, as well as other observables 
such as the strange hadron elliptic flow, will be useful in understanding the dynamics 
of strange quarks and hadrons in heavy ion collisions at RHIC. 

2. A MultiPhase Transport Model - AMPT 

The AMPT model is a hybrid model that uses the minijet partons from hard processes 
and excited strings from soft processes in the HIJING model ^HI f° r the initial 
condition of relativistic heavy ion collisions. The time evolution of partons is then 
described by the ZPC JT] parton cascade model, and that of hadrons by an extended 
ART model ^21- I n the default AMPT model, only minijet partons are included 
in the parton cascade with a parton scattering cross section of a p = 3 mb. After 
partons freeze out, they combine with their parent strings and then fragment to 
hadrons according to the Lund string fragmentation as implemented in PYTHIA • 
The default model has been quite reasonable in describing the measured rapidity 
distributions of charge particles |14l 115) . particle to antiparticle ratios 15 , and the 
spectra of low transverse momentum pions, kaons multi-strange baryons |17j and 
(j> mesons ^S] in heavy ion collisions at SPS and/or RHIC. 

Since the initial energy density in Au+Au collisions at RHIC is expected to be 
much larger than the critical energy density at which the transition from hadronic 
matter to quark-gluon plasma would occur |191 I2(J| . the AMPT model has been 
extended to allow the conversion of initial excited strings to partons at RHIC energies 
|21l [§]. In this string melting scenario, hadrons that would have been produced 
from string fragmentation are converted instead to valence quarks and/or antiquarks. 
Interactions among these partons are again described by the ZPC parton cascade 
model. The transition from the partonic matter to hadronic matter at parton freeze- 
out is achieved using a simple quark coalescence model by combining two nearest 
partons into mesons and three nearest partons into baryons (or anti-baryons) |21|. 
With the energy in excited strings taking part in the early partonic interactions and 
using quark coalescence to model hadronization, the extended AMPT model with 
string melting |22 is able to describe the observed elliptic flow at RHIC [23 I23| . 
which the default AMPT model failed to reproduce. 

At present, the ZPC parton cascade JT] in the AMPT model includes only 
two-parton elastic scatterings. The in-medium differential cross section is given by 
d<jp/di = 9na 2 (l + fj, 2 /s)/2/(t — /i 2 ) 2 , where the effective screening mass \x in principle 
depends on the temperature and density of the partonic matter [TT]. In this study, we 
take \x as a parameter in order to study the effect of partonic scatterings. Also, for 
simplicity, we assume the same scattering cross section for partons of different flavors. 
We note, however, that comparisons of high-quality data on the elliptic flow of strange 
hadrons [25] with theoretical predictions (211 E3 1281 1291 |3"U] is expected to provide 
very useful information on the interactions of strange quarks in dense matter. 

3. Two Ways of Extracting Radius Parameters 

To evaluate the two-kaon correlation function requires the knowledge of the single 
kaon emission function S(x, p). In the AMPT model, it is obtained from the kaon 
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space-time coordinate x and momentum p at kinetic freeze-out. The HBT correlation 
function for two identical hadrons of momenta pi and P2 is then given by ^ 131) 

r CO K\ = 1 J d 4 x 1 d i x 2 S(x 1 ,K)S(x 2 ,K)cos[Q ■ (x 1 - x 2 )} 
2l ^' ' fd*x 1 S(x 1 ,p 1 )Jd*x 2 S(x 2 ,p 2 ) ' {1) 

where K = (pi + p 2 )/2 and Q = (pi — p 2 ,Ei — E 2 ). Expecting that the emission 
function is sufficiently smooth in momentum space, one can evaluate the correlation 
function by using pi and P2 for K in the numerator of above equation. 

First, the size of the emission source can be extracted from the emission function 
via the curvature of the correlation function at Q = 0: 

R%{K) = D Xi , Xj (K) - D Xiif3jt (K) - D 0it ^{K) + D Mt {K). (2) 

In the above, (3 — 'K./Kq with K$ denotes the average energy of the two kaons; 
Xi(i = 1 — 3) are spatial coordinates of a kaon at freeze-out; and D XjV — (x ■ y) — (x) (y) 
with (x) denoting the average value of x. In this study we use the usual "out-side-long" 
(osl) coordinate system j^UHl- 

On the other hand, the measured correlation function C2(Q,K) is usually fitted 
by a Gaussian function in Q, i.e., 



C 2 (Q,K) = 1 + Acxp (^-J2RUK)Q^ , 



(3) 



We note that, for central heavy-ion collisions, the above fitted radius parameters would 
be identical to the source radii given by the curvature of the emission function at Q = 
(i.e. Eq. @) only when the emission source is Gaussian in space and time. 

4. AMPT Results at RHIC Energies 

4--1- Two-kaon correlation function 
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Figure 1. Correlation functions for K° (-K y < 1, 200 < p T < 400 MeV/c). 

Using the program Correlation After Burner |32j . we have evaluated the 
correlation function C2(Q,K) of two kaons in their longitudinally comoving frame 
for central (b — fm) Au+Au collisions at y^s = 130 A GeV. In Figure 1, we 
show the invariant correlation function, i.e., the correlation function as a function 
of Qinv = \Z~Q 2 , and its projections onto one of the Qout, Qsidc and Qi on g axes. In 
evaluating the projected correlation function, the other two Q-components have been 
integrated over the range — 40 MeV/c. The dash-dotted curves in Figure 1 represent 
results from the default AMPT model (no string melting), while the other curves are 
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from the extended AMPT model with string melting but using different values for o~ p . 
It is seen that the one-dimensional kaon correlation functions become narrower as a p 
is increased in the extended AMPT model. Their dependence on a p seems, however, 
to be much weaker than that of the pion correlation function [§]. 

4-2. Source radii versus fitted radius parameters 

The source radii for kaons within — 1 < y < 1 and 200 < pt < 400 MeV/c are shown 
in Figure 2 (a). The results are obtained from both the default AMPT model (shown 
at a p = 0) and the extended AMPT model with string melting. For the latter, we 
have used different parton cross sections of a p = 3, 6, 10, and 15 mb. It is seen that 
these radii have values between 2 and 5 fm, and they are much smaller than the source 
radii for low px pions (between 7 and 25 fm) 9 . The i? ut /-Rside ratio (solid curves 
without symbols) from the kaon emission function is also shown in Figure 2 (a), and 
it increases appreciably with increasing parton scattering cross section. 




Figure 2. (a) Source radii, (b) fitted radius parameters and A for as functions 
of Op at 130A GeV. Points at a p = correspond to the default AMPT. 

The radius parameters obtained from fitting the three-dimensional correlation 
function C2(Q) by Eq. © are shown in Figure 2 (b). We find that they are close 
(mostly within 30%) to the source radii determined directly from the emission function. 
This is contrary to the pion case, where the pion source radii can be more than twice 
larger than the radius parameters from the Gaussian fit UJ. Note that kaons from cf> 
meson decays have been included in the evaluation of the correlation function C*2(Q) 
but not in the source radii (due to the long lifetime of </>); however, pions from u 
decays have been included in the source radii in Ref . [HJ ■ Figure 2 (b) also shows that 
-Rout/-Rside from fitted radius parameters changes little with parton cross section, in 
contrast to the significant increase seen in the source radii shown in Figure 2 (a). 

Figure 3 (a) shows the extended AMPT results with a p — 10 mb for kaon source 
radii as functions of m T in the six p T -bins, - 200 - 400 - 600 - 800 - 1000 - 1500 
GeV/c, for Au+Au collisions at y/s = 130A GeV. We see that for all mr the source 
radii (solid) are close to the corresponding fitted radius parameters (filled squares). We 
note that, before making quantitative predictions on kaon correlations, abundances 
of strange resonances such as K* and 4> need to be checked, and the effects of other 
resonances not yet included in the AMPT model also need to be investigated. 
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Figure 3. (a) Source radii and fitted radius and A parameters for Kg with 
— 1 < y < 1 as functions of mx, as well as fitted radius and A parameters for 
mid-rapidity pions !> . (b) x ou t — iCside and (c) x ut — t distributions at freeze-out 
for Kg. The curve with open diamonds represents {x ou t)(t). 



4-3. x — t correlation and the -Rout/^sidc ratio 

Figure 3 (b) and (c) show the x ou t — ^side and x ou t — t distributions at freeze-out for 
Kg with 200 < pt < 400 MeV/c and — 1 < y < 1. The kaon emission source shows 
a positive shift in x out as in the pion case UJ, consistent with a strong transverse 
flow. The emission source also appears to be highly non-Gaussian, which leads to 
different values of radius parameters extracted from the two methods in Section [3] 
|3*51 • The solid curve with open diamonds in Figure 3 (c) gives the average value 
(x ou t) as a function of freeze-out time t. The kaon emission source is seen to have a 
strong positive a; ou t — t correlation as in the pion emission source 
Since i? s 2 ido = D Xsido , Xsido , but 

-Rout = D x ont ,x out - 2 D Xaut j 3±t + Dp ±tj/3±t , (4) 
the ratio i? ou t / i? S idc contains information about the duration of emission (in the last 
term) and has been studied extensively @J |^ |SJ OS]. We note, however, that a 
direct relation between i? out and the emission duration exists only if the a; ou t — t 
correlation term D Xaat ^ ±t is small, which is not the case according to our results 
from the extended AMPT. E.g, the above equation for the kaons is numerically 
written as 3.4 2 ~ 35 — 2 x 22 + 20 in units of fm 2 . The x ut — t correlation term 
Ac ou t,/3j_t is thus positive (+22 fm 2 ) and comparable to the magnitude of Dp ± t,p ± t 
(20 fm 2 ), making it difficult to extract information about the duration of emission 
from i? ut/-Rsidc- The situation is similar for the pion emission source, e.g., for mid- 
rapidity pions with 125 < px < 225 MeV/c, numerically the corresponding Eq. Q is 
17 2 ~ 185 - 2 x 168 + 431 in units of fm 2 . 

4-4- Energy dependence 

Figure 4 shows the correlation functions for both neutral and charged kaons within 
-1 < y < 1 and 200 < p T < 400 MeV/c from AMPT model with string melting and 
<Tp = 10 mb for \fs = 13QA and 200A GeV. The effect due to Coulomb interactions 
is included for K + correlation functions using the program Correlation After Burner 
[52] . The kaon correlation functions are found to change only slightly from 130A to 
200,4 GeV at RHIC. 
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